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Abstract 

We develop a theory of fluctuations for Brownian systems with weak long-range in- 
teractions. For these systems, there exists a critical point separating a homogeneous 
phase from an inhomogeneous phase. Starting from the stochastic Smoluchowski equa- 
tion governing the evolution of the fluctuating density field of the Brownian particles, we 
determine the expression of the correlation function of the density fluctuations around a 
spatially homogeneous equilibrium distribution. In the stable regime, we find that the 
temporal correlation function of the Fourier components of the density fluctuations de- 
cays exponentially rapidly with the same rate as the one characterizing the damping of a 
perturbation governed by the deterministic mean field Smoluchowski equation (without 
noise). On the other hand, the amplitude of the spatial correlation function in Fourier 
space diverges at the critical point T = T C (or at the instability threshold k = k m ) imply- 
ing that the mean field approximation breaks down close to the critical point and that the 
phase transition from the homogeneous phase to the inhomogeneous phase occurs sooner. 
By contrast, the correlations of the velocity fluctuations remain finite at the critical point 
(or at the instability threshold). We give explicit examples for the Brownian Mean Field 
(BMF) model and for Brownian particles interacting via the gravitational potential and 
via the attractive Yukawa potential. We also introduce a stochastic model of chemotaxis 
for bacterial populations generalizing the deterministic mean field Keller-Segel model by 
taking into account fluctuations and memory effects. 



1 Introduction 

In a recent series of papers [TJ [2], [3l HJ E] , we have considered some theoretical aspects of the 
dynamics and thermodynamics of systems with weak long-range interactions [6], [7] . In these 
systems, the interaction potential u(r) decays with a rate slower than l/r d at large distances, 
where d is the dimension of space (these potentials are sometimes called "non-integrable" ) . 
As a result, any particle feels a potential dominated by interactions with far away particles 
(i.e. the interaction is not restricted to nearest neighbours) and the energy is non-additive. 
This can lead to striking properties (absent in systems with short-range interactions) such as 
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inequivalence of statistical ensembles and negative specific heats in the microcanonical ensem- 
ble. On the other hand, the usual thermodynamic limit N — > +00 with N/V fixed is not 
relevant for these systems and must be reconsidered. If we write the potential of interaction as 
u(|r — r'|) = ku(\r — r'|) where k is the coupling constant, then the appropriate thermodynamic 
limit for weak long-range interactions corresponds to N — > +00 in such a way that the coupling 
constant k ~ 1/N — > and the volume V ~ 1. In that limit, we have an extensive scaling of 
the energy E ~ N and entropy S ~ N (while the temperature T ~ 1), but the system remains 
fundamentally non-additive. Other equivalent combinations of the parameters are possible to 
define the thermodynamic limit as discussed in [H [5] and in various contributions of [7J . For 
systems with weak long-range interactions, it is often claimed that the mean field approxima- 
tion is a very good approximation and that it becomes exact in the proper thermodynamic 
limit N — > +00. In fact, this is true only if we are far from a critical point. Close to a critical 
point, the fluctuations become large and cannot be ignored. In that case, the two-body cor- 
relation function does not factor out in a product of two one-body distribution functions and 
the mean field approximation breaks down. It is therefore highly desirable to derive stochastic 
kinetic equations that go beyond the mean field approximation and that take full account of 
fluctuations. This is the main object of the present paper. 

In our previous studies, we have distinguished two types of systems: Hamiltonian and Brow- 
nian. Hamiltonian systems with long-range interactions are isolated and evolve at fixed energy. 
The dynamics of the particles is described by N coupled deterministic Newton equations. Since 
the energy is conserved, the relevant statistical ensemble is the microcanonical ensemble. The 
evolution of the iV-body distribution function is governed by the Liouville equation and the 
statistical equilibrium state is described by the microcanonical distribution. Examples of such 
systems are provided by stellar systems [HI El UTil HH [12] , two-dimensional vortices [T3 | [T ^ [T5 | [16] 
and the Hamiltonian Mean Field (HMF) model [T71 [18] ; see also the important list of refer- 
ences in these papers. Brownian systems with long-range interactions, on the other hand, are 
dissipative and evolve at fixed temperature. The particles are subject to their mutual long- 
range interactions but they experience, in addition, a friction force and a stochastic force which 
mimic the interaction with a thermal bath (that is due to short-range interactions). Therefore, 
the dynamics of the particles is described by N coupled stochastic Langevin equations. The 
temperature is defined through the Einstein relation as the ratio between the diffusion coeffi- 
cient and the friction coefficient, and it measures the strength of the stochastic force. Since the 
temperature is fixed, the relevant statistical ensemble is the canonical ensemble. The evolution 
of the iV-body distribution function is governed by the Fokker-Planck equation and the sta- 
tistical equilibrium state is described by the canonical distribution. Examples of such systems 
are provided by self-gravitating Brownian particles [121 ED], bacterial populations experiencing 
chemotaxis [HI El 1231 EI] and the Brownian Mean Field (BMF) model [TgllSo] . 

Systems with long-range interactions have a very peculiar dynamics and thermodynamics 
[HI [TJ. When the interaction is attractive, there exists a critical point separating a spatially 
homogeneous (gaseous) phase from a spatially inhomogeneous (clustered) phase. In the micro- 
canonical ensemble (MCE), the "clustered" phase appears below a critical energy E c and in the 
canonical ensemble (CE), it appears below a critical temperature T c . The homogeneous phase 
may still exist for E < E c or T < T c but is dynamically and thermodynamically unstable. 
One fundamental illustration of this type of phase transitions corresponds to the famous Jeans 
instability in astrophysics [261 M- For the Jeans problem 0, the critical temperature T c = +00 

'■As is well-known, the standard Jeans analysis assumes that an infinite and homogeneous self-gravitating 
system is a steady state of the hydrodynamical equations (the so-called Euler-Poisson system), which is not 
correct. This inconsistency in the stability analysis is referred to as the Jeans swindle [9j. However, the Jeans 
treatment can be justified in cosmology |27j if we take into account the expansion of the universe because this 
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so that the system is always in the "clustered" phase: this corresponds to the universe that 
we know, filled of galaxies. A spatially homogeneous distribution of matter is always unsta- 
ble to sufficiently large wavelengths, above the Jeans length kj 1 = (AnGnm 2 /3) -1 / 2 , and this 
leads to gravitational collapse and clustering. If we formally consider a screened Newtonian 
interaction (attractive Yukawa potential) with screening length k^ 1 , we find that the phase 
transition occurs at a finite critical temperature T c = AwGnm 2 / k 2 , separating a homogeneous 
phase from a clustered phase pQ. For T > T c the homogeneous phase is stable and for T < T c 
it becomes unstable to wavenumbers k < k m = ko(T c /T — l) 1 / 2 . Alternatively, we can consider 
a spatially inhomogeneous self-gravitating gas in a spherical box of radius R. In that case, 
there exists a critical energy E c = —0.335GM 2 /R in MCE (discovered by Antonov [281 ES] for 
stellar systems) and a critical temperature T c = GMmj {2.b2ksR) in CE (discovered by Emden 
[30] for isothermal stars and recently emphasized by the author for molecular clouds in contact 
with a thermal bath [31] and for self-gravitating Brownian particles [19]) below which the sys- 
tem passes from a slightly inhomogeneous gaseous phase to a highly inhomogeneous condensed 
phase H. Interestingly, the analogue of these "gravitational phase transitions" also takes place 
in the context of the chemotaxis of bacterial populations in biology [32"l [33] . In these systems, 
an infinite and uniform distribution of cells is a steady state of the equations of motion (the so- 
called Keller-Segel model and its generalizations) so there is no "Jeans swindle". Furthermore, 
the screened Yukawa potential enters naturally in the problem and the screening length has a 
clear physical interpretation as it takes into account the degradation of the secreted chemical. 
Another illustration of this type of phase transitions is given by the study of toy models like 
the Hamiltonian Mean Field (HMF) and the Brownian Mean Field (BMF) models [13 HB]. 
For these systems, there exists a critical energy E c = kM 2 /(8ir) and a critical temperature 
T c = kMj (47r) separating a spatially homogeneous (non-magnetized) phase from a spatially 
inhomogeneous (magnetized) phase. The magnetization plays the role of the order parameter 
and the phase transition is second order in that case [T71 ITS] . We expect the above-mentioned 
types of phase transitions, and their generalizations [3U [33, [3U [371 ESI EE21 EE], to occur for a 
large class of mean field systems with different potentials of interaction. 

The critical point separating the homogeneous phase from the inhomogeneous phase is often 
called a spinodal point. The instability threshold of the homogeneous phase can be obtained 
from different classical methods: (i) by studying the sign of the second order variations of the 
thermodynamical potential (entropy in MCE or free energy in CE) and determining when the 
spatially uniform distribution becomes an unstable saddle point (see, e.g., Sec. 4.4. of Paper I 
or Appendix C of [24]) (ii) by studying the bifurcation (from homogeneous to inhomogeneous) 
of the solutions of the meanfield integrodifferential equation (1-19) determining the statistical 
equilibrium state (see, e.g., Sec. 2.3 of Paper I and Appendix C of [24]) (iii) by studying 
the linear and nonlinear dynamical stability of the homogeneous state with respect to kinetic 
equations: Vlasov, Euler, Landau, Kramers, Smoluchowski... (see, e.g., Paper II and [331 [24] ). 
In the homogeneous phase, the one-body distribution function is trivial (it is spatially uniform 
with a Maxwellian velocity distribution) and the state of the system is usually characterized by 
the two-body distribution function. For weak long-range interactions, the two-body correlation 
function h(\ri — r 2 |) can be obtained from the equilibrium BBGKY-hierarchy at the order 
0(1/N) by neglecting the three-body correlation function that is of order 0(1/N 2 ) [T]. This 
is similar to the Debye-Huckel approximation in plasma physics. It is found however that, for 

creates a sort of "neutralizing background" in the comoving frame expanding with the system allowing for 
infinite and homogeneous steady states. 

2 In order to have a well-defined condensed phase, one has to introduce a small-scale regularization of the 
gravitational potential. More fundamentally, we can invoke quantum mechanics (Pauli exclusion principle) and 
consider the case of self-gravitating fermions (e.g. white dwarfs, neutron stars and fcrmion balls) |12) . 
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large but fixed N, the Fourier transform of the correlation function diverges at the critical 
point (or at the instability threshold) so that the mean field approximation breaks down close 
to the critical point. This implies that the phase transition should take place strictly before the 
critical point (or strictly before the instability threshold) predicted by mean field theory. In 
[H [TS], we have reached this conclusion from the study of the equilibrium BBGKY hierarchy. 
In the present work, we would like to complement our previous studies by developing a theory 
of fluctuations starting directly from the stochastic kinetic equation governing the evolution 
of the fluctuating density field. In this paper, we restrict ourselves to the case of Brownian 
systems. 

The paper is organized as follows. In Sec. [2J we consider a gas of Brownian particles in 
interaction in an overdamped limit where inertial effects are neglected. In Sec. 12.11 using a 
BBGKY-like hierarchy, we give the deterministic kinetic equation (jSJ) satisfied by the smooth 
density profile and the Smoluchowski equation (|Sj) resulting from a mean field approximation. 
In Sec. 12.21 we give the stochastic kinetic equation (ITT)) satisfied by the exact density dis- 
tribution (expressed in terms of 5-functions) and, averaging over the noise, we recover the 
equation obtained from the BBGKY-like hierarchy. In Sec. 12.31 we give the stochastic kinetic 
equation (|2"Tj) satisfied by the coarse-grained density distribution obtained by averaging the 
exact density distribution over a small spatio-temporal window (thus keeping track of fluctu- 
ations). In Appendix [B], we provide another derivation of this equation by using the general 
theory of fluctuations exposed by Landau & Lifshitz [39]. In Sec. 12.4) we derive generalized 
stochastic Cahn-Hilliard equations in the limit of short-range interactions. In Sec. 12.51 starting 
from the stochastic Smoluchowski equation (|44p - (j45|) . we develop a theory of fluctuations for 
Brownian particles with weak long-range interactions. Specifically, we study the correlations 
(<5 Pk(t)5 pk> (t + t)) of the Fourier transform of the density fluctuations and show that it behaves 
like A(k)e a{k)r . In the stable regime, the decay rate cr(k) < of the temporal correlations 
coincides with the decay rate of the perturbations 5pk(t) ~ e CT( - fc ^ governed by the determin- 
istic mean field Smoluchowski equation (without noise). On the other hand, for T < T c , the 
amplitude of the correlation function A(k) diverges as we approach the instability threshold 
k — > k m , suggesting that the instability takes place sooner than predicted by mean field theory 
(for T > T c , considering the "dangerous" mode k = fc*, the correlation function diverges as 
we approach the critical point T — ► T 6 + ). These results clearly demonstrate that fluctuations 
cannot be ignored close to a critical point. In Sec. El we generalize our results to the case of 
an inertial model of particles in interaction including a friction force and a stochastic force. 
The overdamped model is recovered in a strong friction limit. Again, the correlations of the 
density fluctuations diverge as we approach the instability threshold but, in sharp contrast, the 
correlations of the velocity fluctuations remain finite at this threshold. In Sec. 13.41 we derive a 
generalized Smoluchowski equation taking into account memory effects and make contact with 
the Cattaneo model and the telegraph equation. In Sec. 13.51 we derive a stochastic model 
of chemotaxis generalizing the deterministic mean field Keller-Segel model by keeping track of 
fluctuations. Finally, in Sec. HI we extend our study in phase space taking full account of the 
inertia of the particles. We derive a stochastic Kramers equation and make the connection with 
the diffusive and hydrodynamic models of the previous sections. 
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2 The overdamped case 



2.1 The smooth density distribution 



We consider an overdamped system of N Brownian particles in interaction whose dynamics is 
governed by the coupled stochastic equations (see Paper II): 



dt 



(1) 



where /i = l/(m£) is the mobility (£ denotes the friction coefficient), U (ri, rjy) = J2i<j u i\ r i~ 
Tj\) is the potential of interaction, is the diffusion coefficient and Rj(t) is a white noise such 
that (Ri(t)) = and (R?(t)R?(t')) = 5ij5 a p5{t - t 1 ) where i = 1,...,N refer to the particles 
and a = l,...,d to the coordinates of space. We assume that the particles interact via a 
binary potential u(\ti — Tj\) depending only on the absolute distance between the particles. 
As discussed in the Introduction, this stochastic model can describe self-gravitating Brownian 
particles [El [20], bacterial populations experiencing chemotaxis [221 (231 [23] (see also Sec. 13.51) 
or toy models like the BMF model [18]. The iV-body distribution Pjy(ri, rjy, t) is solution 
of the Fokker-Planck equation 

dP N ^ d 3P 



dt 



d_ 



i=l 



N 



<9r ?; 



fxm Pntt- ■ 



(2) 



The stationary solution of this equation is the Gibbs canonical distribution 

1 



P 



N 



-/3m 2 U 



(3) 



provided that the inverse temperature (3 = 1/(/cbT) is related to the mobility and the diffusion 
coefficient through the Einstein relation 



From the Fokker-Planck equation 
first equation of the hierarchy is 



dPi 
~~dt 



_d_ 
dvi 



" = t- (4) 

we can construct a BBGKY-like hierarchy [2] . The exact 
dPx 



+ fim 2 (N - 1) 



du 



12 



PydTn 



The two-body distribution function can be written as 

P 2 (n, r 2 , t) = Pi(n, t)Pi(r 2 , t) + P 2 '( ri , r 2 , t), 



(5) 



(6) 



where P^ri, r 2 , t) is the correlation function. Let us consider a weak long-range potential of 
interaction ii(r 12 ) = ku(r 12 ) in a proper thermodynamic limit iV — > +oo in such a way that 
k ~ 1/N and V ~ 1. Using scaling arguments, it can be shown that P' 2 = 0(1/N) except close 
to a critical point (see [H [5] for details). Therefore, if we are far from a critical point and if iV 
is sufficiently large, we can make the mean field approximation P 2 (ri,r 2 ,t) ~ Pi(i"i, t)Pi(r 2 , t) 
and we obtain 



dP l 



_d_ 
dr 1 



BP 

D,—^ + / im 2 NP 1 (r,t) 
or i 



du 



12 



dri 



P 1 (r 2 ,t)dr 2 



(7) 
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Introducing the smooth density distribution p(r, t) = NmPi(r, t), this equation can be rewritten 

as 



dp 
di 



[r, t) = D*Ap(r, t) + pmV ■ (p{r, t)V J p(r', t)u(r - r')dr^j 



It can be put in the form of a mean field Smoluchowski equation 

dp 

= A,Ap + pmV- (pV$), 
where <3>(r, t) is a smooth potential produced by the particles themselves 



$(r,t) = J p(r',t)u(r-r') dr'. 



(8) 



(9) 



(10) 



If we introduce the mean field Boltzmann free energy functional 

F 



E-TS=- [ p$ dr + k B T [ L l n L dr, 

2 mm 



we can rewrite the mean field Smoluchowski equation (Q in the form 

V- 



dp 
di 



1 / v. __ SF 



(12) 



This equation satisfies an if -theorem appropriate to the canonical ensemble 



6£dp 
Jpdi 



dr 



6£ 
Jp~ 



V 



1 ^5F 

i op 



dr 



(13) 



The steady solution of the mean field Smoluchowski equation (jH]) or (1121) corresponds to a 
uniform p = 5F/5p leading to the mean field Boltzmann distribution 



p = Ae 



(14) 



where $(r) is given by Eq. (jiop . Finally, we note that the mean field Smoluchowski equation 
(JBJ) can be written in Fourier space as 



^(k, t) = -LU 2 p(k, t) - (2ir) d {im [ k • k'p(k - k', t)u(k')p(k', t)dk'. 
dt J 



(15) 



For weak long-range potentials of interaction, the mean field approximation usually provides 
a good and useful description of the system as a first approach □. We must, however, recall 
its domain of validity: (i) first, it assumes that the number of particles in the system is large 
(mathematically speaking it is valid in the limit iV — > +oo). Therefore, we can expect deviations 
from mean field theory due to finite N effects. These deviations will become manifest for 



3 It is often advocated that long-range potentials of interaction exhibit lack of temperedness and stability. 
Some potentials, like the cosine potential in the HMF model are well-behaved. By contrast, the gravitational 
potential is singular and, strictly speaking, there is no statistical equilibrium state (no global maximum of 
entropy in MCE and no global minimum of free energy in CE). However, there exist long-lived metastable 
states (local maxima of entropy or local minima of free energy) that can be adequately described by the mean 
field approximation [12j . The formation of a Dirac peak in CE, which can be viewed as the "equilibrium state" 
of the system, can also be described by the mean field Smoluchowski-Poisson system [40]. By contrast, the 
formation of binary stars in MCE requires going beyond the mean field approximation. 
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sufficiently large times (see discussion at the end of Sec. I2.3p . (ii) Close to a critical point 
T —> T c , the correlation function diverges (see pQ and Sees. 12.51 and [2761) . Typically, we expect 
a scaling of the form P' 2 ~ N~ X {T — T c )~ l so that the limits N — > +00 and T — > T c do not 
commute (see Sec. 2.7 in [IB]). Therefore, even for large N, the mean field approximation 
is expected to break down close to a critical point because P' 2 is not necessarily small (the 
mean field approximation is valid for iV ^> (T — T c ) _1 , which requires larger and larger particle 
numbers as T — > T c ). In the two cases (i) and (ii) mentioned above, we must take fluctuations 
into account and consider stochastic kinetic equations as discussed in the sequel. 



2.2 The exact density distribution 



The exact density distribution of the particles is expressed as a sum of Dirac distributions in 
the form 



N 



[r,t) = m^25(r-ri(t)). 



(16) 



i=l 



It was shown by Dean [H] (see also P21B31IH] ) that the exact density field satisfies a stochastic 
equation of the form 



dp d 
dt 



[r,t) = D*Ap d (r,t) + pmV ■ U*(r,t)V J p d (r',t)u(r - r')dA 

+V- (y2D m mp d (r,t)B.(T,t) 



(17) 



where R(r, t) is a Gaussian random field such that (R(r, £)) = and (R a (r, t)R^(r', t')) = 
SapS(r — r')5(t — t'). Note that the noise is multiplicative [H]. Introducing the exact potential 



$ d (r,t) = J p d (r',t)u(r-r') dr', 



the stochastic equation (fT7|) can be rewritten as 

^ = D*A Pd + /imV • (p d V$ d ) + V ■ (y/2D*m Pd R). 

For = and p 7^ 0, we get the exact deterministic equation 

dpd 



dt 



pmV ■ (pdV$d). 



If we introduce the discrete Boltzmann free energy functional 

1 



E d — TSa 



p d §d dr + k B T 



Pd-. Pd , 
— In — dr, 

m m 



(19) 



(20) 



(21) 



we can rewrite the stochastic equation ffl9l in the form 



dpd 
dt 



V 



1 ( +\x? 5Fd 
-p d {r,t)V— 



V 



2k B Tp a 



R 



(22) 



This equation can be viewed as a Langevin equation for the field p d (r,t). From this equation, 
it can be shown that the probability of the density distribution t] is governed by the 

Fokker-Planck equation 



dW[p d ,t] 
dt 



S 



Spd(r,t) 



V- p d (r,t)V 



D*— + p— 

opd op d 



W[p d ,t}\dr. 



(23) 
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At equilibrium, we get the Boltzmann distribution W[p d ] oc e -P(Fd[Pd]-»f Pddr) glj S3J, IB] . 

If we average Eq. (JTT1) over the noise, we find that the evolution of the smooth density field 
p(r, t) = (p d ) is governed by an equation of the form 

^(r, t) = £*Ap(r, t) + fimV ■ J (p d (r, t)p d {r\ i))Vu(r - r')dr' . (24) 

Using the identity (see Appendix lAl) : 

(p d (r, t)p d (r', t)) = Nm 2 P 1 (r 1 t)5(r - r') + N(N - l)m 2 P 2 (r, r', t), (25) 

and assuming that Vtt(O) = 0, we find that Eq. (121"]) coincides with the exact equation ([5]) of the 
BBGKY-like hierarchy giving the evolution of the one-body distribution function. Furthermore, 
if we make the mean field approximation (p d (r, t)p d (r', t)) ~ p(r, t)p(r', t), we recover the mean 
field Smoluchowski equation (jHJ). 



2.3 The coarse-grained density distribution 

Equation (jSJ) (or equivalently Eq. (124|) ) for the ensemble averaged density field p(r,t) is a 
deterministic equation since we have averaged over the noise. In contrast, Eq. (1171) for the 
exact density field p d (r,t) is a stochastic equation taking into account fluctuations. However, 
it is not very useful in practice since the field p d (r,t) is a sum of Dirac distributions, not a 
regular function. Therefore, it is easier to directly solve the stochastic equations ([1]) rather than 
the equivalent Eq. (fT7|) . Following [H], we can keep track of fluctuations while avoiding the 
problem of 5-functions by defining a "coarse-grained" density field p(r, t) obtained by averaging 
the exact density field on a spatio-temporal window of finite resolution. The "coarse-grained" 
density field satisfies a stochastic equation of the form 



dp 
dt 



>, t) = D*Ap(r, t) + fimV ■ (J p (2) (r, r, i)Vu(r - r')dr'j 



+V • (y/2D*mp(T,t)R(r, t)) . (26) 

where ~p^ (r, r', t) is a two-body correlation function. For a weak long-range potential of interac- 
tion and for a sufficiently small spatio-temporal window, we propose to make the approximation 
p( 2 )(r, r', £) ~ p(r, t)p(r', t). In that case, we obtain a stochastic equation of the form 

^(r, t) = D*Ap(r, t) + pmV ■ (p(r, t) J p(r', t)Vu{r - r')dr'J 

+V • ( v / 2 J D,mp(r,t)R(r, t)) . (27) 

Introducing the smooth potential 

$(r, t) = J p(r, t)u{r - r) rfr', (28) 
the stochastic equation 02 7p can be rewritten 

^ = D,Ap + pmV • (7jV$) + V ■ (y/2D*mp R). (29) 

This will be called the stochastic Smoluchowski equation for the smoothed-out density field 
p(r, t). This equation is intermediate between Eqs. (jHJ) and (117j) . It keeps track of fluctuations 
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while dealing with a continuous density field instead of a sum of 5-functions. This equation 
will be central in the rest of the paper. We will see that it can reproduce the equilibrium 
density correlation function (1701) that was obtained in Paper I from the equilibrium BBGKY- 
like hierarchy p]. Therefore, it represents an improvement with respect to the mean field 
Smoluchowski equation flS]). We stress that this equation is physically distinct from Eq. (TTTj) . 
In Appendix |Bj we propose an alternative derivation of Eq. (1291) by using the general theory 
of fluctuations exposed by Landau & Lifshitz |39j . 

If we introduce the coarse-grained Boltzmann free energy functional 



This equation can be viewed as a Langevin equation for the field p(r, t). The evolution of 
the probability of the density distribution W[p, t] is governed by a Fokker-Planck equation of 
the form ( 1231) where Fa and pd are replaced by F cg , and p. At equilibrium, we have W[p] oc 
e -/3(F c . g .[p]-fifpdr)^ p or jy- _^ .(.q^ ^he equilibrium distribution W[p] is strongly peaked around 
the global minimum of F c _ g , [p] at fixed mass. However, the system can remain trapped in a 
metastable state (local minimum of F c , g \p\) for a very long time. Let us be more precise. If 
we ignore the noise term, Eq. ( 13TT) reduces to Eq. ( 1121) . In that case, the system tends to a 
steady state that is a minimum (global or local) of the free energy functional _F c . ff . [p] at fixed 
mass (maxima or saddle points of free energy are linearly dynamically unstable with respect 
to mean field Fokker-Planck equations [21]). If the free energy admits several local minima, 
the selection of the steady state will depend on a notion of basin of attraction. Without noise, 
the system remains on a minimum of free energy forever. Now, in the presence of noise, the 
fluctuations can induce a dynamical phase transition from one minimum to the other. Thus, 
accounting correctly for fluctuations is very important when there exists metastable states. The 
probability of transition scales as e _AF / fcsT where AF is the barrier of free energy between two 
minima. Therefore, in an infinite time, the system will explore all the minima and will spend 
most time in the global minimum for which AF is the largest. Now, for systems with long-range 
interactions, the barrier of free energy AF scales as N so that the probability of escape from 
a local minimum is very small and behaves as e~ N . Therefore, even if the global minimum is 
in principle the most probable state, metastable states are highly robust in practice since their 
lifetime scales like e . They are thus fully relevant for N ^> 1 [17]. These interesting features 
(basin of attraction, dynamical phase transitions, metastability,...) would be interesting to 
study in more detail in the case of systems with long-range interactions. Some results in this 
direction have been reported in [15J SSI HI] in the gravitational case. 

2.4 Generalized Cahn-Hilliard equations 

Let us assume that w(|r — r'|) is a short-range potential of interaction and that the preceding 
equation fl2Ul) remains valid (to simplify the notations, we drop the bar on the coarse-grained 
fields). Setting q = r' — r and writing 




(30) 



we can write the stochastic equation ( 1291) in the form 




(31) 




(32) 
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we can Taylor expand p(r + q, t) to second order in q so that 

, . . . dp 1 v-^ d 2 p 

p(r + q, t) = p(r, t) + ^ — 9i + - ^ ^fcfc. 



Substituting this expansion in Eq. (1321) . we obtain [ISj: 

$(r,t) = -ap(r,t) - ^Ap(r,t), 



(33) 



(34) 



with a = —Sd Jq°° u(q)q d ~ 1 dq and b = — ^Sd J +oc u{q)q d+l dq. Note that / = {b/a) 1 / 2 has the 
dimension of a length corresponding to the range of the interaction. Substituting Eq. 
Eq. (l30l . we can put the free energy in the form 



in 



FcM = J 
where V is the effective potential 

V{p) = - 

In that case, Eq. ( 13~TT) can be rewritten 
dp 



-(Vp) 2 + V(p) 



a 2 2k B T 

-p H — plnp + Vq. 

b mb 



^ = -AV-\pV(Ap-V'(p))]+V- 



2k B Tp 



R 



with A = 6/(2£). Substituting Eq. (13"6T) in Eq. (13T1) we explicitly obtain 

= ^Ap - ^Ap 2 - b -V ■ (pV(Ap)) + V • (y^Tfr R 



(35) 



(36) 



(37) 



(38) 



Without the noise term, the steady state of Eq. (l3Tj) . (I37|) or (1381 corresponds to a uniform 
p = 5F/ 5p leading to 

2a 2k B T , 2k B T , , 

Ap = - p = --p + -5- lnp + -5- - p. (39) 



In particular, at T = 0, Eq. (1381) reduces to 



f| = -|A/-|v.(pV(Ap)), 
and its steady state is solution of the Helmholtz equation 



2a 



-P. 



(40) 



(41) 



Morphologically, Eq. (l3Tj) with Eq. (1351) . or equivalently Eq. (1371) . is similar to the stochastic 
Cahn-Hilliard equation [49J for model B (conserved dynamics): 



e^ = A^ + v / 2|^TV-R, F[p] = J 



-(Vp) 2 + V{p) 



or explicitly 



dp 
dt 



-A(Ap - V'(p)) + ^2£k B TV ■ R. 



(42) 



(43) 



There are, however, crucial differences between Eqs. ( 13?]) and (143]) . First, the presence of the 
density p(r, t) in the deterministic current and in the noise term. Secondly, in the usual Cahn- 
Hilliard problem, the potential has a double-well shape of the typical form V(p) = (1 — p 2 ) 2 
leading to a phase separation while, in the present case, the potential (I3"B"|) is of a different 
nature. 
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2.5 Theory of fluctuations 



Let us return to the stochastic Smoluchowski equation (131]) satisfied by the coarse-grained 
density distribution that we write in the form (for convenience, we drop the bars on the coarse- 
grained fields): 



dt 



1 (k B T 



£ V m 



- Vp + pV$ 



j 2k B Tp 



R 



(44) 



$(r,t) = J p(r',i)u(r - r) dr' 



(45) 



We wish to study the fluctuations of the density around an infinite and homogeneous equilibrium 
distribution. To that purpose, we consider small perturbations Sp(r,t) and 5<3>(r, t) around the 
steady state p(r) = p, $(r) = $ with $ = p f u(x)<ix. The linearized equations for the 
perturbations are 



d8j) = kBT A5 + M$ + ^bJ^N ■ R, 
at m 



(46) 



5$(r, t) = y Sp(r' ,t)u(r — r') dr' . 

We now decompose the perturbations in Fourier modes in the form 

5p(r,t) = J Spfcuy^-^dkdu, 

and similar expressions for 5$(r, t) and R(r, t). Taking the Fourier transform of Eqs. 
f|4"7|) . we obtain the algebraic equations 



i^u5p 



kui 



rn 



k 2 5p ku) - pk 2 5$ ku) + ^2k B TipV ■ R 



ku)i 



(47) 



(48) 



61) and 



(49) 



5$ kuJ = (27r) d u(k)8p kwi 



(50) 



where we have used the fact that the integral in Eq. (1471) is a convolution. Solving for Sp ku} , we 
obtain 



k B T , 



k 2 + (2n) d u(k)k 2 p - i& 



Sf>ku 



2k B Tip ik^R 



where the correlations of the Fourier transform of the noise are given by 



(51) 



(52) 



Without noise (R = 0), Eq. ( IBTT) gives the dispersion relation (see Paper II) associated with 
the mean field Smoluchowski equation ([9]), i.e.: 



/ T 1 

Z(k, u) = —k 2 + (27r) d u(k)k 2 p - i£u = 0. 
m 



(53) 
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The function Z(k,u) plays a role similar to the dielectric function in plasma physics. The 
perturbation evolves exponentially rapidly like Sp k (t) oc e "^* with a rate given by a(k) = 
—ujQ^k)/^ where 

u; 2 (k) = -^—k 2 + (2n) d u(k)k 2 p. (54) 
m 

Thus, we find that the spatially homogeneous phase is stable with respect to the mean field 
Smoluchowski equation (jUJ) if u 2 {k) > 0' i- e - 

1 + (2n) d ppmu(k) > 0, (55) 

for all modes k and unstable (to some modes) otherwise. If u > 0, the homogeneous phase is 
always stable. Otherwise, a necessary condition of instability is that 

k B T < k B T c = {2ix) d pm\u{k)\ max . (56) 

If this condition is fulfilled, the range of unstable wavenumbers is determined by 

7 rp 

(2ir) d \u(k)\ > —. (57) 
pm 

Some explicit examples of potentials of interaction, and the corresponding conditions of insta- 
bility, are given in Paper I and in [33J. 

Let us now determine the correlations of the fluctuations around a stable equilibrium ho- 
mogeneous distribution in the presence of noise. If we take the noise into account (R^ 0), the 
Fourier transform of the density fluctuations is given by 



Sp k . = — ^ . (58) 



yJ2k B T^p iWR? 
Z(k, uj 

Therefore, the correlations of the fluctuations in Fourier space are 



ma sh ) - -^BT^pk^(RlR^) 



Using Eq. (1521) . we get 



(Sp^SM = ^yj^p i* + 10% + a/), (60) 



or, more explicitly, 



(Spku$Pk'u') = \rr\ TT^ 2k BT£,pk — ^ + k')5(a; + uj'). (61) 



(27r) d + l \k^L k 2 + {2n) d u{k)k 2 p\ 2 + i 



2 oo 2 



The temporal correlation function of the Fourier components of the density fluctuations is given 
by 

{5p k (t)5p k ,(t + t))= [ (6p k Jp k ^)e wt e w '^ +T) dujduj'. (62) 



Using Eq. (IBTI) . the integral on to' is trivial and yields 

(5p k (t)6pk>(t + r)> = —^2k B T^6{\s. + k') J du (63) 
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The integral on ui can be performed by using the Cauchy residue theorem. The poles of the 
integrand are the zeros of the functions Z(k,ou) and Z(k,u)*, i.e. they are solutions of the 
dispersion relation (153]) and its complex conjugate. If u%(k) > 0, which is required for the 
stability of the homogeneous phase, the integrand has a single pole in the upper-half plane at 
to = 2c<Jq(/c)/£ and the residue is e~ u °( k > T >t /[2i£ujQ(k)]. Therefore, we obtain 

(Sp k (t)Sp k '(t + r)) = (a^^f ^ + kOe- 1 ^^, (64) 
or, more explicitly, 

(Spk{t)5p k ,{t + r)) = 5{k + k')e L« p \ /? . (65) 

( 27r ) + {2n) d u(k)p 

This formula is one of the most important results of this paper. We shall come back to its 

physical interpretation in Sec. I2.61 The equal time correlation function (corresponding to 
r = 0) is given by 

»' W ^ uWH t ) i|t + 4 (66) 

In the absence of any interaction between the particles (u = 0), we get 

(SpkSpk') = j^ypm5(k + k'), (67) 
which corresponds to the standard result [39J: 

<(M 2 > = ff (68) 
On the other hand, using Eq. ( l66l) and the identity (see Appendix lAl) : 

(SpkSpki) = J^yP™ I 1 + W l nh(k)j 5(k + k'), (69) 
we find that the Fourier transform of the correlation function is 

eq[ } ~ l + (27r) d Pnm 2 u(k)- [ ' 

This is precisely the result (1-54) obtained in Paper I by analysing the second equation of the 
equilibrium BBGKY-like hierarchy and neglecting the three-body correlation function. There- 
fore, the stochastic Smoluchowski equation (I44|) - (l45|) is able to reproduce the equilibrium two- 
body correlation function. On the other hand, from Eqs. (j65l) and (J69l) . the Fourier transform 
of the equilibrium temporal correlation function is 

h(k,t,t + r) = ^(k)e-[^ fe2+(2w) ^ )fc2 ^ /? . (71) 

This can be compared to the out-of-equilibrium temporal evolution of the equal-time spatial 
correlation function h(k, t) given by Eq. (11-165) of Paper II. Note that the condition of stability 
(155]) is implied by Eq. (163]) . Eq. fTTTj) . Eq. (11-165) and Eq. (155"]) . This completes the discussion 
given in Sec. 4.2 of Paper I. 
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2.6 Specific examples 

The physical content of formula (164)) is very instructive. First, we note that the temporal cor- 
relation function of the Fourier components of the density fluctuations decreases exponentially 
rapidly with a decay rate cr(fc) = —u^k)/^ that coincides with the decay rate of a perturbation 
of the density governed by the deterministic mean field Smoluchowski equation ([§]), i.e. without 
noise 0. According to this temporal factor, or according to the mean field theory based on the 
Smoluchowski equation ([9]), the modes satisfying the inequality ui^k) > should be stable. 
For T < T c , the threshold of instability corresponds to the wavenumber(s) k = k m where k m is 
defined by u:l{k m ) = 0. Now, we note that the amplitude of the fluctuations in formula (1541) 
behaves like o>q(/c) _1 , so that it diverges as we approach the instability threshold k — > k m . On 
the other hand, if we denote by k* m the value of the critical wavenumber at T = T c satisfying 
1 + (2ir) d ftcpmuik^n) = and if we consider the mode k = k* m in Eq. (I65I) . we get for T > T c : 

mt)Sp»{t + r)> = ^-lEL 6{k + k> -MT-T c)(fc * ) V (^) (72) 

This formula clearly shows that the correlation function diverges at the critical point T = T c for 
the "dangerous" mode k = k* m . These results imply that the mean field approximation breaks 
down close to the critical point (or close to the instability threshold) and that the instability 
triggering the phase transition can occur sooner than what is predicted by mean field theory 
(i.e. from the stability analysis of the mean field Smoluchowski equation). Some results in this 
direction have been reported in [451 S3 HZ] in the gravitational case. 

Let us consider specific examples for illustration (we use the notations of Paper I). For 
the BMF model (the Brownian version of the HMF model) [18] . there exists a critical tem- 
perature T c = kM/(4n). Considering the linear dynamical stability of a spatially homoge- 
neous distribution with respect to the mean field Smoluchowski equation (11961) . we find that 
ujQ(n) = Tn 2 + 2iru n pn 2 where u n is given by Eq. (I198p . The modes n^il decay exponentially 
rapidly as e~ Tn t ^ so they are always stable. By contrast, the modes n = ±1 evolve in time 
like e~ <yT ~ Tc ^ t ^ . For T > T c , the perturbation is damped and for T < T c the perturbation has 
exponential growth. In that case, the homogeneous phase is unstable to the n = ±1 modes (see 
Appendix [C]) . According to formula (164"|) . the correlation function of the density fluctuations 
can be written for the stable modes n ^ ±1: 

(5p n (t)5p m (t + r)> = J^- e - Tn2 ^5 n ,. m , (73) 
and for the "dangerous modes" n — ±1: 

M T 

(Sp ±1 (t)Sp m (t + r)> = ^J—Yf iT ~ Tc)T/ ' S rn, T i- (74) 

This simple toy model, where the potential of interaction is restricted to one Fourier mode, 
is very interesting for pedagogical purposes because it clearly illustrates the discussion given 
above. Considering the temporal factor in Eq. (1741) . we see that the correlations decay for T > 
T c with the rate given by mean field theory. However, as we approach the critical temperature 
from above (T — > T c + ), the amplitude of the fluctuations diverges like (T — T c )~ l implying that 

4 Similarly, in Sec. 2.9 of Paper II, we noted that, for Hamiltonian systems with long-range interactions, 
the Fourier modes of the temporal correlation function of the force decay exponentially rapidly with a decay 
rate that coincides with the decay rate of a perturbation of the distribution function governed by the Vlasov 
equation. 
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a ( k) = = ^_JL^[^ {Te/T _ i) _ (75) 



the mean field approximation breaks down [f| and that the phase transition should occur for T 
strictly above T c . We had also reached this conclusion in [18] from the study of the equilibrium 
BBGKY-like hierarchy. 

Let us now consider the attractive Yukawa potential pp. A detailed stability analysis of 
the homogeneous phase with respect to the mean field Smoluchowski equation (and general- 
izations) has been performed in [33] • There exists a critical temperature ksT c = S^Gprajk^ 
depending on the screening length k$ l . Furthermore, in the linear regime, the perturbation 
evolves exponentially rapidly as Spk(t) oc e a( - k ^ t with a rate 

u 2 (k) _ k B T k 2 
£ m£ k 2 + &o 

For T > T c , the homogeneous phase is always stable and for T < T c it is unstable to wavenum- 
bers k < k m (T) = ko(T c /T — l) 1 / 2 . Taking into account the fluctuations, formula (IMl) shows 
that the correlation function of the density fluctuations is 

(5 Pk (t)5 Pk , (t + T)) = j^y k2 + k , {1 _ TJT) *(k + k )e ^ . (76) 

Considering the mode fc = fc^ = 0, we see that the correlation function diverges like (1— T c /T) _1 
as we approach the critical temperature T — > T c + . On the other hand, for T < T c , we see that 
the amplitude diverges like {k 2 — k"^)" 1 as we approach the critical wavenumber k — > fc+(T). 
This is particularly true for the gravitational interaction (ko = 0) for which Eq. (!76|) reduces 
to 

mt)8p»<t+T)) = ( 2^p- fc ^ ( k+k/ ) e "^ (fc2 " fc2j)f ' ( 7? ) 

where fcj = (SdGmPp) 1 ^ 2 is the Jeans wavenumber. According to the standard Jeans analysis 
[26], the homogeneous phase is stable against perturbations with wavenumbers k > kj and it 
becomes unstable for k < kj. However, the divergence of the correlation function as k — > k^j 
suggests that the gravitational instability will take place sooner, i.e. for smaller wavelengths 
than the Jeans length. This conclusion was previously reached by Monaghan [15] on the basis 
of a hydrodynamical model of self-gravitating system incorporating viscosity and fluctuations. 



3 The inertial model 

In this section, we generalise the previous results to the case of a stochastic model taking into 
account inertial effects. This corresponds to the damped Euler equations with a long-range 
potential of interaction and a stochastic forcing. This generalization allows us to study the 
correlations of the fluctuations of the velocity field and their behaviour close to the critical 
point. The stochastic Smoluchowski equation (144j) - (145j) is recovered in a strong friction limit 
£ — > +oo by neglecting the inertial term (l.h.s.) in Eq. (1751) . The stochastic damped Euler 
equations can find applications in certain biological models of chemotaxis where inertial effects 
are relevant [50], EH [23] (see also Sec. 13. 5p . 

3.1 The density correlations 

We consider the stochastic damped Euler equations 

^ + V ■ (pu) = 0, (78) 
5 This implies that the limits N — > +oo (mean field) and T — > T c do not commute. 
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<9u / „s 
^ + (u-V)u 



k B T 



m 



Vp - pV$ - £pu - ^2k B T£pIl(r, t), 



(79) 



$(r,t) = J u(r-r')p(r',t) dr', 



(80) 



where the quantities have their usual meaning. Without the stochastic term (R = 0), we recover 
the damped Euler equations introduced in [52] (see also [201 [231 ESI El] ) • With the stochastic 
term, we obtain a more general model taking into account fluctuations. The form of the noise is 
justified in Appendix [B] using the general theory of fluctuations of Landau & Lifshitz [39]. For 
£ = 0, we get the usual Euler equations and for £ — > +oo, neglecting the inertial term (l.h.s.) in 
Eq. (1791 and substituting the resulting expression £pu ~ — {k B T '/ 'm)V p — pV '$ — ^2k B T^pH in 
the continuity equation (1781) . we recover the stochastic Smoluchowski equation (l4^ - (l4"5l) . Note 
that Eq. ( 1791 can be written in terms of the free energy fTTTT) as 



P 



du 

~dt 



[u- V)u 



,SF 



-pV-^ - £pu - v/2^T^R(r, t). 
5p 



Considering small perturbations around a uniform distribution with p(r) 
and u = 0, we find that the linearized equations for the perturbations are 



d5p 
~dt 



pV • u = 0, 



p, $(r) = $ 
(82) 



k B T. 



u = _w_ v8 _ pW$ _ ^ _ ^ 2 fc B TfpR(r, t), 
at m 



(83) 



<J$(r, t) = y u (r - r')Sp(r', t) dr'. 



They can be combined to give 



c/t 2 at m 



34) 



(85) 



If we decompose the perturbations in Fourier modes of the form e^ k r ut >, we obtain the system 
of algebraic equations 



iuSpkuj + ipk • u ku> = 0, 



(86) 



lUpVL ku 



k B T 



m 



k5p kuJ - ikp5$ ku - £pu fcw - ^/2k B T£pR 



kuii 



»7) 



<5$ feaJ = (27r) d n(A;)5p^. 
Solving for the density perturbation, we get 



knT 



m 



k 2 + (2n) d u(k)pk 2 - uo(u + i£) 



Spk. = y/2k B T£p iWR^. 



(88) 



(89) 
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Without noise (R = 0), Eq. (|89p gives the dispersion relation associated with the mean field 
damped Euler equation [20], i.e. 

7 T 1 

Z(k, u) = — k 2 + (2n) d u(k)pk 2 - lu(uj + if) = 0. (90) 
m 

Using the definition (154|) . the dispersion relation can be rewritten 

u 2 + i£u - u 2 {k) = 0, (91) 



with solutions 

^H=^f±M«,-4 ±!W . (92) 

The perturbation evolves in time as e - " 1 ^^ : . liuj 2 (k) < 0, the perturbation grows exponentially 
rapidly in time with a rate [— £ + a/£ 2 + 4|o^q |]/2. If uj 2 (k) > f 2 /4, the perturbation decays 
exponentially rapidly in time with a rate — f/2 and oscillates with a pulsation ^/ 4cJq — f 2 /2. 
If < ^o(^) < C 2 /4, the perturbation decays exponentially rapidly in time with a rate [— f + 
a/£ 2 — 4cg>q]/2 without oscillating. Therefore, according to mean field theory, the system is 
stable iff oo 2 (k) > for all k. This returns the condition fl55l) studied in Sec. 12.51 For T < T c , 
the onset of instability corresponds to the wavenumber(s) k m such that Uo(k m ) = 0. 

We now consider stable modes {oj 2 {k) > 0) and study the correlations of the density fluc- 
tuations in the presence of noise. Repeating the steps of Sec. 12.51 going from Eq. (155]) to Eq. 
(|63|) . the density correlation function can be written 

(6p k (t)6p k ,(t + T)) = ^-1^2k B npk 2 I(k,r)S(k + k'), (93) 

where I(k, r) is the integral 

r+oo i 

f(cv)e^du, f(u) = (94) 

\Z(k,uj)\ 2 

We can evaluate the integral with the Cauchy residue theorem. The poles of f{w) are the zeros 
of the functions Z(k,u) and Z(k,u)*, i.e. they are solutions of the dispersion relation fl9~Tj) and 
its complex conjugate. If uj 2 (k) > 0, the function f(uj) has two simple poles in the upper-half 
plane at u = if/2 + Q and uj = if/2 — Q. The residues of f(oj)e lu)T at uj = if/2 + Q and at 
uj = if /2 — Q are 



if (if + 2ft) (2fi) ' (if - 2ft)zf (-2ft) 
Using the Cauchy residue theorem, and recalling that 



ft 2 = c 2 - if 2 , (96) 



we obtain after simplification 



1 = T^ e 



f 

cos(ftr) + — sin(ftr) 



(97) 
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where Q can be complex. In conclusion, the density correlation function for the inertial model 
is 



(Sp k (t)Spk'(t + r)) 



1 k B Tpk 2 _ in 
-e 2 



(2kY w g(fc) 



cos(f2r) 



2Q 



sin(f2r) 



<5(k + k') 



(98) 



The interpretation of this formula is similar to the one given in Sec. 12.61 In particular, we see 
that the amplitude of the fluctuations diverges when k — > k m , corresponding to uj 2 (k) — > 0. 
Thus, instability should set in slightly before we reach the range of unstable wavenumbers 
(157]) . In particular, the amplitude diverges at the critical point T — > T c + for the "dangerous" 
wavenumber k* m . For £ — »■ (Euler), we can make the approximation Q ~ cj , and we obtain 



(8pk(t)6p k .(t + T)) 



1 A; B TpA; 2 _£ 



e 2 T cos(u;oT)^(k + k'). 



For £ 



n 



(2nY u 2 (k) 

+oo (Smoluchowski), we can make the approximations 



-i 1 



2^ 2 



e 2 



cos(f2r 



l _^d. 

~ -e « 
2 



1 2 

_£ J- "o . 

e 2 r sin(fir) ~ e £ ' 

v ; 2i 



2" v" e 

and we recover the result of Eq. (j64|) obtained in the overdamped limit. 



(99) 



(100) 



3.2 The velocity correlations 

Let us now turn to the correlations of the velocity fluctuations. From Eq. (!H7|) . the Fourier 
components of the velocity fluctuations satisfy the relation 



p(£ - iu})u kuJ = --j^ik5p kL0 - \/2k B Tip R kiAj . (101) 
Therefore, the velocity correlations can be expressed as 

p 2 {i - - vJ) (utA) = ~j^ ktlk ' V QPkJM + 2k B n P (RlRl, wl ) 

2 '2 

+^ik^2k B Tip{5p kuj Rl, wl ) + ^ik' v ^2k B Ttp{5p k ^Ri u ). (102) 



Using the relation 



iVWT P k a R a kul 

Opku = TTTj x , (103) 

Z(k, LO) 



and Eq. (J52|) . we get 



^4 



hi 2 k^k a ~ - ii) 2 k v k a 

^"T-zm^^ ~ 2kBT ^W {ieM (104) 



or, more explicitly, 



p 2 (e + u 2 ) (ulu^,) = -^2k B T£p6(k + k')S(u; + u/) 

x ™, _ ^o 2 _ "1 WW \ 

\ k 2 \Z(k,uj)\ 2 k 2 Z(k,u) k 2 Z{k,tu)*j' { ' 
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Finally, using the identity 



1 1 2Re(Z) 

V/L , .\* 7/1. . . \ I 2 IV/L ,.M2 1 l iUD J 



Z(Jfe,w) Z(A;,cu)* \Z(k,uj)\ 2 \Z(k,u)\ 
the foregoing relation can be rewritten 

fc 2 |Z(£;,^)| 2 P |Z(fc,cu)| 2 J ' 1 ' 

Taking the inverse Fourier transform in cu-space of this relation, we find that the temporal 
correlations of the velocity fluctuations are given by 

{mm + r)> = {2n)d+1 B p M k + k') {-^f-^K + ^f-^K' + JW J , (108) 

where we have introduced the integrals 

„ V*MP' A " = L V*MP' (109) 



A"' 



do;- 



110) 



These integrals are easily calculated with the Cauchy residue theorem. After simplification, we 
obtain 



1 k R T 



(27r) d p 



A; 2 cu 2 + 2£ 2 



1 



— 2 — 1 ) cos(f2r) 







KH) sin(nT) -( 3+ 5H1 +e M' (in) 



We note that, at variance with the density correlation function 
function does not diverge when k —>■ k m . Indeed, for Uq = 0, we have 



1 k B T 
(27r) d p 



e^ T 5(k + k')^. 



the velocity correlation 



;ii2) 



On the other hand, taking r = in Eq. (11 111) , we obtain the equal time velocity correlation 
function 



1 k B T 



(2vr) d p 



-5(k + k') 5" u 



Auj 2 k»k y 
out + 2£ 2 k 2 



Contracting the indices, we get 

(U fc (t) • Ufc'(t)) 



(2ir) d p 

For £ — > (Euler), Eq. (jllip can be simplified into 

1 k B T 



1 kBT 5(k + k') 2de ~ {A ~ d)uJl 



(27T) d p 



-5(k + k'] 



k^k u _i 

p 2 ' 

k 2 

19 



co 2 + 2e 



cos(a;or) + 3e 2 r 



113) 



114) 



115) 



Alternatively, for £ — > +00 (Smoluchowski), we get 

«(tK,(t + r)) = - (27r)d I flk + lfl^— e (116) 

This result can be obtained directly from the study of the stochastic Smoluchowski equation in 
Sec. 12.51 by defining the velocity field 

7 T 1 

£pu = — — Vp - pV$ - y/2k B T£pR. (117) 

m 

Note that the velocity correlations in Eq. fl 1 1 6 [) tend to zero when k — > k m contrary to Eq. 
(11121) . This shows that the limits £ — > +00 and k — > fc m do not commute. 

3.3 Specific examples 

Let us discuss specific examples by restricting ourselves, for brevity, to the overdamped limit 
(11161) . For the stochastic BMF model, we get for the stable modes il: 

QT2 2 

(u n (t)u m (t + r)) = —^^-ne- TnTli , (118) 
and for the "dangerous" modes n — ±1: 

ST 

+ ^)> = --^^,ti(T " T c )e-( T " T ^. (119) 

The velocity correlations tend to zero when T — > T 6 + . For the attractive Yukawa potential, we 

get 

(120) 

For T < T c , the amplitude tends to zero as we approach the critical wavenumber k — > fc+(T). 
For the gravitational interaction (ko = 0), Eq. (11201) reduces to 

+ 0> = - (i^(* 2 - *5)^e-^" ! -5)f i(k + 10. (121) 

The amplitude goes to zero as k —>■ kj. The fact that the velocity correlation function does 
not diverge when k — > kj was previously observed by Monaghan [15] with his hydrodynamic 
model. 

3.4 Stochastic model with memory 

The stochastic damped Euler equations (1781) - (1791 can be rewritten 

^ + V-(pu) = 0, (122) 

|-(pu) + V(pu ® u) = Vp - pV$ - £pu - y/2k B T^pK(r, t). (123) 

at m 
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If we neglect the inertial term (l.h.s.) in Eq. (11231) and substitute the resulting expression for 
pu in Eq. fl 1 22j) . we obtain the stochastic Smoluchowski equation (JHJ). This is valid in a strong 
friction limit £ — > +00. We can obtain a more general model taking into account some memory 
effects. If we neglect only the nonlinear term V(pu Cg> u) in Eq. (1 1 23j) . we get 

d , s kuT 



-(pu) = — ^-Vp - pV$ - ipu - y^ffpRir, t). (124) 

Taking the time derivative of Eq. (11221) and substituting Eq. (I124p in the resulting expression, 
we obtain a simplified stochastic model keeping track of memory effects 



d_P f dp 
8t 2 5 dt 



V • (^Vp + pV^j + V • [^2k B T^pR) . (125) 
In terms of the free energy (fill) , we have 

^4r^<) +v '(^w)- 

This equation, which is second order in time, is analogous to the telegraph equation which 
generalizes the diffusion equation by introducing memory effects. We note that in the linear 
regime |u| <C 1 considered in Sec. El the nonlinear term V(pu(g>u) in Eq. ( 11231) is negligible so 
that Eq. (11251) can be justified rigorously from the damped Euler equations in this regime. This 
implies that the theory of fluctuations that we have developed in Sec. |3] directly applies to the 
stochastic equation (11251) . In particular, the linearization of Eq. (11251) around a homogeneous 
distribution returns Eq. (1851) . However, the stochastic Smoluchowski equation with memory 
(11251) may also be relevant in the nonlinear regime as a heuristic equation. Indeed, although we 
have neglected the nonlinear term V(pu ® u) in Eq. (I123p . we have kept the full nonlinearities 
in the right hand side. Therefore, Eq. (11251) is a semi-linear model intermediate between the 
fully nonlinear hydrodynamical model (1 1 2 21) - (1 1 23 j) and the linearized hydrodynamical model 

(USD. 

Finally, we note that Eq. (11261) can be viewed as a form of stochastic Cattaneo model. The 
deterministic Smoluchowski equation can be written as a continuity equation dtp = — V- J where 
the current J = — (l/£)pVp is proportional to the gradient of a chemical potential p = 5F/5p 
[24J. This is similar to Fick's law for the diffusion of particles or to Fourier's law for the 
diffusion of heat. In the context of heat conduction, Cattaneo [53] has proposed a modification 
of Fourier's law in order to describe heat conduction with finite speed. He assumed that the 
current is not instantaneously equal to the gradient Vp but relaxes to it with a time constant 
1/t. In the present situation, this would lead to a model of the form 

^ + V-J = 0, (127) 
83 1 5F hk B Tp 

where we have included the stochastic term for completeness. These equations are equivalent 
to the semi-linear model formed by Eqs. (11221) and ( 11241) if we set J = pu and r = l/£. For 
t = 0, we recover the stochastic Smoluchowski equation (13~T1) . More generally, taking the time 
derivative of Eq. (1 1 2 T[) and using Eq. (1 1 28j) . we obtain 

^ + t-V.(ipV^Vv.|,/^R|, (129) 



R, (128) 



dt 2 dt V.r 5p J 

which coincides with Eq. ( 11261) provided that we take r 
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3.5 Application to chemotaxis 



In this section, we briefly mention the application of the preceding results to the problem of 
chemotaxis in biology |54j . A more detailed discussion is given in a specific paper [55] with 
complements and amplification. The standard Keller-Segel (KS) model [21] of chemotaxis can 
be viewed as a form of mean field Smoluchowski equation [24]. It describes the diffusion of 
bacteria (or other chemotactic species) in the concentration gradient of a chemical produced by 
the particles themselves. As we have seen in this paper, the correlation function diverges close 
to a critical point. In that case, the mean field approximation breaks down and the fluctuations 
must be taken into account. Fluctuations also play an important role when the particle number 
./V is small and when there exist metastable states (local minima of free energy). In that case, 
fluctuations can trigger dynamical phase transitions from one state to the other (see Sec. 12. 3j) . 
For these different reasons, it is important to derive a chemotactic model going beyond the 
mean field approximation and taking into account fluctuations. 

We start from a microscopic model of chemotaxis where the dynamics of the particles is 
governed by iV coupled stochastic equations of the form 

^ = xVc d ( ri (t), t) + y/2D~*Ri(t), (130) 



?^ = -kc d + D c Ac d + h^26(T-T i (t)), (131) 

i=l 

where r^t), with i = 1, N, denote the positions of the particles and c<j(r, t) is the exact field 
of secreted chemical. In these equations, x an d -D* represent the mobility and the diffusion 
coefficient of the organisms and k, h and D c represent the degradation rate, the production 
rate and the diffusion coefficient of the secreted chemical. By extending Dean's approach (and 
the results of Sec. 12.31) to the case of chemotactic species, we obtain a stochastic Keller-Segel 
model of chemotaxis: 



dp 



r, t) = £>,Ap(r, i) - X V ■ (p(r, t)Vc(r, t)) + V ■ (s/2D,p(r,t)R(r, t)) , (132) 

dc 

— (r, t) = -kc(r, t) + D c Ac(r, t) + hp(r, t), (133) 

generalizing the deterministic mean field Keller-Segel model. This model fully takes into account 
the effect of fluctuations 0. On the other hand, there exists situations in biology where inertial 
effects must be taken into account [50]. In that case, parabolic models like the Keller-Segel 
model must be replaced by hyperbolic models similar to hydrodynamic equations [501 EH E3] . 
By extending the results of Sec. [3] we obtain a hydrodynamic model of chemotaxis taking into 
account inertial effects and fluctuations in the form: 

H + V ■ {pa) = 0, (134) 



6 Note added: Until now, fluctuations have been ignored by people working on chemotaxis. Therefore, our 
paper is the first attempt to include fluctuations in the Keller-Segel model. However, after submission of this 
paper arXiv:0803.0263], a paper by Tailleur & Cates [arXiv:0803.1069] came out on a related subject. These 
authors also consider the effect of fluctuations in the motion of bacteria. They derive transport coefficients 
from microscopic models but do not take into account the long-range interaction between bacteria due to 
chemotaxis. Alternatively, in our approach, the coefficients -D* and x appearing in the Langcvin equations 
arc phcnomcnological coefficients but chemotaxis is fully taken into account. Therefore, these two independent 
studies are complementary to each other. 
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— {pa) + V(pu <g> u) = -f D.Vp + pVc - £pu - v/2D,f 2 pR(r, i), (135) 

coupled to the field equation (11331) . In the strong friction limit £ — > +00 where the inertial 
term in Eq. H135[) can be neglected, it returns the stochastic KS model H132[) with x = !/£• On 
the other hand, if we only neglect the term V(pu ® u) in Eq. (I135P like in Sec. 13.41 we obtain 
a stochastic equation of the form 

X^? + % = V ■ {D*Vp - xpVc) + V ■ (v/^R) , (136) 

It can be viewed as a stochastic Cattaneo model of chemotaxis (or a stochastic telegraph 
equation) . 



4 The stochastic Kramers equation 

In this section, we generalize the results of Sees. 12.1112.31 in phase space. This is the rigorous 
way to take into account inertial effects and fluctuations in the problem. The motion of the 
Brownian particles is described by N coupled stochastic Langevin equations of the form (see 
Paper I): 

(137 » 

^ = -£vi - mW l U{v u v N ) + V2DRi(t). (138) 

The friction coefficient £ and the diffusion coefficient D are related to each other by the Einstein 
relation £ = Df3m where (3 = l/(/c#T) is the inverse temperature [I]. In the strong friction 
limit £ — > +00, we can neglect the inertial term in Eq. fl 1 3 8 1) and we obtain the overdamped 
equations (CQ) of Sec. Owith p = l/(fm) and D, = D/£ 2 . 

Extending Dean's approach [JT] in phase space, we find that the exact distribution function 
/d(r, v, t) — m J2iLi <5( r— r i(^))^( v— V «W) expressed in terms of 5-functions satisfies a stochastic 
equation of the form 



(139) 



where Q(r, v, t) is a Gaussian random field such that (Q(r, v, t)) = and (Q«( r ) v ; t)Q/3( r 'i v '; 0) 
= 5 a p5{r — v')8(y — v')<5(t — t') and $d(r, £) is defined by Eq. (fTBl . If we average over the noise, 
we obtain 

f + V ' ^ ~ t ' / rfrW '[ VM ( r " r ')] ^ v ' ^ r '> v ' *)> = t " (^f + e/V ) • (140) 

Using / = NmPi and the identity 

(fd(r, v, t)/ d (r', v', t)) = Nm 2 P 1 (r, v, t)<f(r - r')5(v - v ') + iV(iV - l)m 2 P 2 (r, v, r', v', t),(141) 

we find that Eq. (j!40p is equivalent to Eq. (11-139) obtained from the BBGKY-like hierarchy. 
If we implement a mean field approximation (/^(r, v, i)/<i(r', v', £)) ~ /(r, v, t)f(r', v', i), we 
obtain the mean field Kramers equation [2]: 
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where $(r, t) is defined by Eq. (fTUl) . Finally, we can heuristically propose a stochastic kinetic 
equation for the evolution of the coarse-grained distribution function f(r,v,t) obtained by 
averaging fd{r, v, t) over a small spatio-temporal window. This leads to the stochastic Kramers 
equation 



df df - df d f df \ d 



where $(r, t) is defined by Eq. (1281) . This equation keeps track of fluctuations but applies 
to a continuous distribution function instead of a sum of Dirac distributions. An altnernative 
derivation of this equation is proposed in Appendix [B] using the general theory of fluctuations 
of Landau & Lifshitz [39] . 

Let us now try to make the link with the parabolic and hydrodynamic models considered in 
Sees. |5] and [31 Taking the hydrodynamic moments on the stochastic Kramers equation (11391) 
and proceeding as in [S21 ED] , we obtain 

^ + V ■ (pu) = 0, (144) 

' (t + = -§| - "S - - / » < i45 > 

where p(r,t) = J fdv is the density, u(r, i) = (1/p) / /vdv is the local velocity, w = v — 
u(r, i) is the relative velocity and Pij = J fwiWjdv is the pressure tensor. Defining g(r, t) = 
J ^2DmfQdv, it is clear that g is a Gaussian noise and that its correlation function is 

(g l (r,t)g J (r' 1 t')} = 2Dm J y/f(r, v, t)/(r', v', t')(Q t (r, v, t)Qj{v\ v', t'))dvdv' 

= 2Dm6 ij 6{r - r')5{t - t') J /(r, v, t)dv = 2Dm5 ij 5(r - r')6{t - t')p(r, t). (146) 

Therefore, the equation for the momentum (11451) can be rewritten 

f dm dui\ dPa <9$ . /— „ . r ^ 

{w +u ^) = -^- f '^ i - (pu - ~ ^^•< r - «)■ < 147 > 



p 



This equation is not closed since the pressure tensor depends on the next order moment of 
the velocity. If, following [52j[20], we make a local thermodynamic equilibrium (L.T.E.) ap- 
proximation /iTE( r )V, t) ~ (Pm/27r) d / 2 p(r,t)e~ /3mw I 2 to compute the pressure tensor, we find 
that ~ (k B T/m)p5ij. In that case, Eqs. (11441) and (11471) yield the stochastic damped Euler 
equations (J78 j) -( [79j) . We recall, however, that there is no rigorous justification for this local ther- 
modynamic equilibrium approximation. Therefore, it does not appear possible to rigorously 
derive the damped Euler equations (1751) - (1791 from the Kramers equation (11431) . Alternatively, 
if we consider the strong friction limit £ — > +oo for fixed /3, leading to D = £/(/3m) — * +oo, the 
first term in the r.h.s. of Eq. (TT39|) implies that f(r,v,t) ~ (f3m/2n) d / 2 p(r,t)e- pmv2/2 + O(l/0, 
u = 0(l/£) and P i& = (k B T/m)p5 ij + 0(1/0- To leading order in l/£, Eq. (fT4"7D becomes 



1 fknT 



Vp + pV$ + V^DmpR(r,t)j . (148) 



Inserting Eq. (11481) in the continuity equation (j!44p and defining p = 1 / (£m) and D* = D/£ 2 = 
k B T / (^m), we obtain the stochastic Smoluchowski equation (1291) . This equation can thus be 
derived from Eq. ( 11431) in the limit £ — > +oo. 
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5 Conclusion 



In this paper, we have developed a theory of fluctuations for a system of Brownian particles 
with weak long-range interactions. Starting from the stochastic Smoluchowski equation (]44ft - 
(I45p . justified in Appendix [B] from the Landau & Lifshitz general theory, we have obtained 
a simple formula (1651) for the temporal correlation function of the Fourier components of the 
density fluctuations at equilibrium (for an infinite and homogeneous distribution). This formula 
shows that the correlations decay in time with the same damping rate as the one obtained from 
the study of the normal modes of the deterministic Smoluchowski equation (Q, without noise. 
Furthermore, the amplitude of the correlation function diverges at the critical point T c (or 
at the instability threshold k = k m ) leading to a failure of the mean field approximation in 
that case. As a result, the limits N — > +00 and T — > T c do not commute and the instability 
occurs strictly before the critical point as discussed in [431 \%Q\ |4"T] for gravitational systems. In 
future works, we shall extend this theory of fluctuations to more general models. Indeed, the 
method developed in this paper can be generalized to any type of kinetic equations including 
fluctuations. In particular, the structure of formula ( |63i) where Z{k,uj) is a sort of "dielectric 
function" obtained from the linearized kinetic equation without noise, has a general scope. 



A Correlation functions 

Considering the correlation function of the exact density field (Tl6l) . and introducing the one 
and two-body distributions, we find that 

N 

( Pd (r)p d (r')) = (m 2 ^ 6(r - r,)5(r' - r,)> = (m 2 ^ 6(r - r,)5(r' - r)> 

i,j i=l 

+ (m 2 S(r - r,i)5(r' - Tj)) = iVm 2 Pi(r)5(r - r') + N(N - l)m 2 P 2 (r, r'). (149) 

Denoting by p(r) = (p(r)) = NmP\{r) the equilibrium averaged distribution and introducing 
the fluctuations Sp(r) = p^(r) — p(r), we get 

(5p(v)5p(v')) = (p d (v)p d (v')) - p(v)p(v'). (150) 

Starting from the identity (11491) and introducing the correlation function h(r, r') through the 
defining relation P 2 (r,r') = P 1 (r)P 1 (r')[l + h(r, r') + l/N], we obtain at the order 0(1/N): 

(5p(r)5p(r')) = mp(r)5(r - r') + p(r)p(r')h(r, r'). (151) 

For a spatially homogeneous equilibrium distribution where p(r) = p = mn and h(r, r') = 
h(\r — r'|), the Fourier transform of the density fluctuations is 

(Wpk') = tA*P™ [l + (27r)^(k)] 8(k + k'). (152) 

The equilibrium correlation function can be obtained from the equilibrium BBGKY-like 
hierarchy using a Debye-Hiickel-type of approximation or from field theoretical methods using 
the Landau approximation. Starting from the Gibbs canonical distribution in configuration 
space (1-44), which is the steady state of the iV-body Smoluchowski equation (T5]), we can obtain 
P] the equilibrium BBGKY-like hierarchy (1-45). The first two equations of this hierarchy are 

9Pl -(N-l)Pm 2 f^P 2 dr 2 , (153) 



dri J dri 
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dP2 2)/?- 2 / *r 3 . (154) 



Introducing the decomposition (1-14) in Eq. (I153p . we first obtain 
dPi . . 2 [ du\i , . s o /" 



Sri 



{N-l)/3m 2 / ^^(^^(^^-(iV-l^m 2 / ^P 2 '(r 1; r 2 )dr 2 . (155) 
J OTi J 9ri 

Then, introducing the decomposition (I- 14)- (I- 15) in Eq. (1154ft . and using Eq. (11551) to simplify 

some terms, we get 

f 2 = -^P 1 {v 1 )P x {v 2 ) 9 ^ -/W^( ri) r 2 )^H 
or i or i or i 

+Pm 2 j ^P 1 (r 1 )P 1 (r 2 )P 1 (r 3 )dr 3 -(N-2)(3m 2 J ^H^( ri , r 2 )P 1 (r 3 )rfr 3 
/ ^ P 2(ri,r3)Pi(r 2 )rfr3-(iV-2)/5m 2 | ^J*(r 2) r 3 )Pi(ri)dr 3 

-(TV - 2){3m 2 J ^Ps(rx, r 2) r 3 )dr 3 . (156) 



-/3m 



At the order 1/iV in the thermodynamic limit defined in Paper I where Pi, /3, m, \r\ are 0(1), 
P 2 ', u are 0(1/N) and P 3 are 0(l/iV 2 ), the foregoing equations reduce to [J 

^J. = _(AT- l)/3m 2 / ^P^rOP^r,)^ - iV/3m 2 / ^P£(r lt r 2 )dr 2 , (157) 
OTi J dri y dri 

9P 2 a 2 r, / \ id i \ du ^ , a 2 f du l3 



dri 



-/5m 2 P 1 (r 1 )Pi(r 2 )^ + /5m 2 / ^P^P^P^dr, 
or i J or i 



-NPm 2 j ^-P^r^P^drs-NPm 2 J ^P 2 '(r 2 , r^r^. (158) 

Now, introducing P^(l,2) = P 1 (l)P 1 (2)[/i(l, 2) + i] in Eq. (1T58|) . using Eq. (fToTj) . and neglect- 
ing terms of order 0(1/N 2 ) or smaller, we find that the correlation function satisfies 

h(n, r 2 ) = -/3m 2 u 12 - N(3m 2 J u 13 h(r 2 , r 3 )P 1 (r 3 )rfr 3 , (159) 

where Pi(r) is given by the zeroth order Eqs. (1-20) and (1-21). We emphasize that this relation, 
which was not given in Paper I, is valid for a possibly spatially inhomogeneous equilibrium state. 
For a homogeneous distribution, we recover Eq. (1-51) which can be solved in Fourier space 
yielding Eq. (l70l . 

It is very instructive to recover these results in a different manner using methods of field 
theory [56]. The equilibrium probability of the density distribution governed by the stochastic 
Smoluchowski equation ([29} is W[p) = ± e -K F W-i* I > r ) with Z = JVp e -fH.m-nI p**) (to sim- 
plify the notations, we drop the bars on the coarse-grained fields). To compute the correlation 
function G(r, r') = (5p(r)Sp(r')), it proves convenient to introduce an auxiliary field /i(r) and 
write 

W[p] = L e -mp]-fn(r)pdr), ( 16Q ) 



7 Note that some terms of order 1/N were missing in Paper I because we made the approximation N — 1 ~ N 
everywhere which is incorrect. 
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The equilibrium corresponds to /i(r) = \i. In the Landau (mean field) approximation, the free 
energy F = —ksTlnZ is given by F ~ F[p] — J p(r)pdr where p(r) is the most probable 
distribution of VT[p]. The maximum of W[p] satisfies the condition /i(r) = 5F/5p(r). Using 
the expression (|30|) of the free energy, we obtain 

u(r — r')p(r')dr' H lnp(r). (161) 

m 

At equilibrium, taking /i(r) = p, we recover the mean field Boltzmann distribution (fT4l) . On 
the other hand, taking the functional derivative of Eq. (I16ip with respect to /i(r') and using the 
fundamental identity G(r,r') = kBT5p(r)/5p(r') [56], we find that the equilibrium correlation 
function is solution of 

5(r -r')=/3 [ u(r - r")G(r', v")dr" + — U-G(r, r'). (162) 
J mp[Y) 

Finally, substituting Eq. (11511) in Eq. (11621) and simplifying some terms, we recover Eq. (I159p . 
Noting that the partition function of the iV-body problem can be written 

Z = J e- f)m2u dr l ...dr N = J Vp(r)e s[p] e-^ E[p] = J Vp(v)e- mp \ (163) 

where the sum runs over the macrostates p(r) with mass J p(r)dv = M and e s ^ denotes the 
number of microstates associated with the macrostates p(r), we see the link between the two 
previously exposed methods. 



B Application of the Landau-Lifshitz general theory of 
fluctuations 

In this Appendix, we derive the stochastic Smoluchowski equation (|44 T) -( l45T) by using the general 
theory of fluctuations exposed in Landau & Lifshitz (see [39], Chap. XVII). We write the 
equation for the density in the conservative form 

| = -V-J, (164) 



where J is the current: 



X f k T \ 



(165) 



The first term is the deterministic Smoluchowski current (see, e.g., [2]) and the second term 
is a stochastic term that takes into account fluctuations. The problem at hand consists in 
characterizing the stochastic term g(r,t). In order to use the general theory of fluctuations 
[39J , we divide the fluid volume in small elements AV and take the average of each quantity 
in each element. The continuum limit AV — > will be performed in the final expressions. 
Equations (11641) and (11651) correspond to the equations 



x n 



S y^ j -y ab X b + y a , (166) 



of the general theory [39] provided that we make the identifications x a — » — J a and y a — > g a . 
The X a can be obtained from the expression of the rate of production of entropy. In fact, since 
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we are working in the canonical ensemble, the proper thermodynamical potential is the free 
energy F = E — TS or, equivalently, the Massieu function J = S — E /T which is the Legendre 
transform of the entropy. It can be written explicitly 

J = -k B [ -^ln-^-dr-i- [ p<$>dv. (167) 
J m m 2 1 J 

Taking the time derivative of this expression, using Eq. (I164p . and integrating by parts, we 
obtain the expression 

j= - / ( — Vp + pV® J • 3 dr. (168) 



Tp \ m 

Note that for g = (no noise) we recover the appropriate form of the H-theorem valid in the 
canonical ensemble [2j [20] : 

j= —- ( -?—Vp + P V$ J dr > 0. (169) 



Tpt; \ m 

If we replace the integral in Eq. (I168P by a summation on AV, we obtain 

J = -£^(^ V ^ V$ )- JA ^- ( 17 °) 

According to the general theory [39], we must also have 

j = -k B ^X a x a . (171) 

a 

Comparing Eq. (11701) with the general expression (I17ip . we find that the X a are given by 



ksTp \ m dx a dx c 

It is now easy to find the expression of the coefficients jab that appear in Eq. (11661) . Comparing 
Eqs. ffT65|) . (HMD and ffT72|) . we find that 

7ab = (if a ^ b); ~j aa = ( 173 ) 

Now, the general theory of fluctuations [39] gives 

(y a (ti)y b (t 2 )) = (7afc + IbaWr - t 2 ). (174) 
Therefore, the correlation function of the stochastic field g{r,t) satisfies 

(g a (r,t)gp(r',t')) = (if r^r'), (175) 

(g a (v,t)g,( r ,t')) = 2 -^5 aP 5{t-t'). (176) 
Taking the limit AV — > 0, we can condense the above formulae under the form 

(g a (r, t)g p (v', t')) = 2 -^5 aP 5{t - t')5{v - r'). (177) 
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We thus recover the expression of the stochastic term appearing in Eq. (I44p by a method 
different from Dean [4Tj . 

We can repeat the same arguments for the inertial model (!78|) - (l80|) . We write Eq. (I79j) in 
the form 



P 



- + (u.V)u 



7 T 1 

; -Vp-pV<&-^u-g'(r,t), (178) 



in 



where g'(r, t) is the noise term to be determined. Comparing Eq. (I178P with Eq. (11661) . we 
have the correspondances x a — > —£,pu a — g' a and y a — > —g' a - The Massieu function J = S — E/T 
for the inertial model is 1201: 



J=-k B [ In L dr - — [ p$ dr - — [ pu 2 dr. (179) 
J m m IT J ' 2T J K ' 

Taking the time derivative of this expression and using Eqs. (I78"]) - (l8"D]) . we obtain after some 
elementary calculations (see, e.g., Appendix G of [21]) the expression 

j = \\ u-(£pu + g')dr. (180) 

For g' = (no noise), we recover the appropriate form of the if -theorem valid in the canonical 
ensemble for the mean field damped Euler equation [20] : 

1 

T 



J = - I £pu 2 dr > 0. (181) 



The discrete expression of Eq. (11801) is 

j = ^u-(^n + g , )AK (182) 

Comparing Eq. (11821) with the general expression (I17ip . we find that the X a are given by 

AV 

X a - —u a . (183) 
k B T 

Then, comparing Eqs. (11781) . (I166P and (11831) . we find that 

lab = (if a ^6); 7aa = ^P. (184) 

Finally, using Eq. (I174p . we obtain the correlation function 

(g' a {T,t)g' p {v', t')) = 2k B T&6 aP 6(t - t')6(r - r'), (185) 



which coincides with the expression given in Eq. (1791) . 

Let us finally briefly consider the kinetic model of Sec. HI It can be written in the form 

df df v$ df d j 
dt dr dv dv ' 



with the current 

J = - D I t 
dv 



D(?f + (3mfv)-g;(r,v,t). (187) 
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The corresponding free energy (Massieu function) J = S — E/T is explicitly given by [2J: 

J=-k B [ — Xn — drdw-^- f f—drdv-^- [ p$dr. (188) 
mm T J 2 2T H 



It production rate is 



J= -k B I -i- ( ^ + /3m/v ) • Jdrdv. (189) 



mf \dv 



For g = (no noise), we recover the appropriate form of the i?-theorem valid in the canonical 
ensemble for the mean field Kramers equation [2J, [20] : 

j = k B [ — s + /3mfv) drdv > 0. (190) 



mf \dv 

On the other hand, repeating the general procedure developed previously, we find that 

X a - ~ ( + Pmfva) AV, (191) 



mf \dv a 

lab = (if a^b); laa = ^J-, (192) 

where AV is the elementary volume in phase space. Passing to the limit AV — > 0, this leads 
to a correlation function of the form 

(g a (r,v,t)g p (r', v',t')) = 2Dmf5 aP 5(t - t')5(r - r')S(v - v'), (193) 

which coincides with the expression given in Eq. ( 11391) . 

C Dispersion relation for the inertial BMF model 

In this Appendix, we complement the discussion of Sec. 13. II by studying the dispersion relation 
associated with the damped Euler equations ( 178l - (l80l) without noise (R = 0) for the inertial 
BMF model [H]. These equations can be written as 

=- r |-sf sin( * (195) 

For £ = 0, they reduce to the Euler equations (see Eq. (67) of [15] ) and for £ — > +oo, we can 
neglect the inertial term in Eq. (j!95p and we obtain the mean field Smoluchowski equation (see 
Eq. (234) of [IB]): 
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Figure 1: Evolution of the dangerous modes n = ±1 for the BMF model with £ = (Euler). 
For T > T c , the homogeneous phase is stable and the perturbation oscillates with pulsation uj. 
For T < T c , the homogeneous phase is unstable and the perturbation increases exponentially 
rapidly with a growth rate 7 > 0. 

Considering the linear dynamical stability of a homogeneous distribution with respect to the 
damped Euler equations (j!94p - (jl95p . and decomposing the perturbation in normal modes 5f ~ 
e i{nQ-u)t) ^ dispersion relation fl90l) can be written 

u(uj + if) = Tn 2 + 27ru n pn 2 , (197) 

where 

k 

U n = - — (5n,l + K-l)- (198) 
47T 

Let us first consider the case £ = (Euler). For n ^ ±1, the dispersion relation becomes 
uj 2 = Tn 2 so that the perturbation oscillates with a pulsation uj = y/T\n\. For n — ±1, the 
dispersion relation becomes uj 2 = T — T c where T c = kMj (4ir) is the critical temperature of 
the BMF model [18]. For T > T c , the perturbation oscillates with a pulsation cj = \/T — T c 
and for T < T c , the perturbation grows exponentially with a gowth rate 7 = \/T c — T. In that 
case, the homogeneous phase is unstable (see Fig. [Q. 

Let us now consider the overdamped case £ — > +00 (Smoluchowski). For n 7^ ±1, the 
dispersion relation becomes i^uo = Tn 2 so that the perturbation decays exponentially with a 
rate 7 = — Tn 2 /£. For n — ±1, the dispersion relation becomes ifu; = T — T c . For T > T c , the 
perturbation decays with a decay rate 7 = — (T — T c )/f and for T < T c , the perturbation grows 
exponentially with a rate 7 = (T c — T)/f . In that case, the homogeneous phase is unstable (see 
Fig. E]). 

In the general case, setting a = —iu so that the perturbations behave as 5f ~ e CTt , the 
dispersion relation (j!97p can be rewritten 

a 2 + f a + Tn 2 + 2-iru n pn 2 = 0. (199) 

For n 7^ ±1, it reduces to a 2 + f cr 4- Tn 2 = 0. The solutions of this equation are <r± = 
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12 3 

T 

Figure 2: Evolution of the dangerous modes n = ±1 for the BMF model with £ — > +oo 
(Smoluchowski). For T > T c , the homogeneous phase is stable and the perturbation decays 
with a rate 7 < 0. For T < T c , the homogeneous phase is unstable and the perturbation grows 
with a rate 7 > 0. 




12 3 4 

T 

Figure 3: Evolution of the dangerous modes n = ±1 for the inertial BMF model described by 
the damped Euler equations. A homogeneous distribution is unstable for T < T c and stable 
for T > T c . For T < T c , the perturbation grows with a rate 7 > 0. For T c < T < T*, 
the perturbation decays with a rate 7 < without oscillating. For T > T*, the perturbation 
undergoes damped oscillations with a decay rate £/2 and a pulsation 
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-£ ± y/ A n )/2 where A n = £ 2 — 4Tn 2 . Let us introduce the wavenumber 

e 2 ^ 1/2 



= ( ^ ) • (200) 

For A n < 0, corresponding to n 2 > n 2 , the perturbation presents damped oscillations with a 
pulsation and decay rate 

co = VT(n 2 - nl) 1 ' 2 , 7 = ~t/2. (201) 

For A n > 0, corresponding to n 2 < n 2 (since we have assumed n ^ ±1, this regime is accessible 
iff In* | > 2 i.e. T < £ 2 /16), the perturbation has a pure exponential decay with a damping rate 

7 = -f + v / T(n 2 -n 2 ) 1 / 2 . (202) 

The modes n 7^ ±1 are always stable, whatever the temperature. For the "dangerous" modes 
n — ±1, the dispersion relation becomes cr 2 + £<r + T — T c = 0. The solutions are a± = 
(— £ ± y/A~x)/2 where Ai = £ 2 — 4(T — T c ). Let us introduce the temperature 

T* = T c + ^. (203) 

For Ai < 0, corresponding to T > T*, the perturbation undergoes damped oscillations with a 
pulsation and decay rate 



00 — a/T — T*, 7 = -£/2. (204) 

For < A x < £ 2 , corresponding to T c < T < T*, the perturbation has a pure exponential decay 
with a damping rate 



I + y/T^T. (205) 



7 2 

For Ai > £ 2 , corresponding to T < T c , the perturbation grows exponentially rapidly with a 
growth rate 



7 = -| + VT,-T. (206) 

Therefore, for T < T c , the homogeneous phase is unstable to the modes n = ±1. The growth 
rate is maximum for T = with value 7* = 7 (0) = — £/2 + V^T- The dependence of the growth 
rate, damping rate and pulsation as a function of the temperature are plotted in Fig. [3j This 
figure should be compared with Fig. 1 of [33] obtained for the gravitational potential (in this 
analogy, the critical temperature T c plays the same role as the Jeans wavenumber k 2 ). 
It can be useful to introduce a dimensionless number 

F - e/A - e - ^ (207) 

F ~ ~t7 ~ w P ~ m> (207) 

which measures the efficiency of the friction force (a similar number has been introduced in 
[201 S3] ) . It can be written as F ~ (££d) 2 where t D ~ \ j\fkp is a typical dynamical time (see 
Sec. 2.2. of [5|). Thus, F is the ratio of the dynamical time on the friction time r ~ l/£. In 
terms of this parameter, the temperature (12031) marking the appearance of oscillations can be 
written T* =T C (1 + F). 
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